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Abstract 

In this work, we review and correct the first Born differential cross section for the pro- 
cess of Mott scattering of a Dirac-Volkov electron, namely, the expression (26) derived by 
Szymanowski et al [Physical Review A 56, 3846 (1997)]. In particular, we disagree with the 
expression of (^^^ they obtained and we give the exact coefficients multiplying the various 
Bessel functions appearing in the scattering differential cross section. Comparison of our 
numerical calculations with those of Szymanowski et al. shows qualitative and quantitative 
differences when the incoming total electron energy and the electric field strength are increased 
particularly in the direction of the laser propagation. Such corrections are very important since 
the relativistic electronic dressing of any Dirac-Volkov charged particle gives rise to these coef- 
ficients that multiply the various Bessel functions and the relativistic study of other processes 
(such as excitation, ionisation, etc..) depends strongly of the correctness and reliability of 
the calculations for this process of Mott Scattering in presence of a laser field. Our work has 
been accepted [Y. Attaourti, B. Manaut, Physical Review A 68, 067401 (2003)] but only as 
a comment. In this paper, we give the full details of the calculations as well as the clear 
explanation of the large discrepancies that their results could cause when working in the ultra 
relativistic regime and using a very strong laser field corresponding to an electric field e = 5.89 
in atomic units. 

PACS number(s): 34.80.Qb, 12.20.Ds 

1 Introduction 

In a pioneering and very often cited paper , Szymanowski et al. [1] have studied the Mott scattering 
process in a strong laser field. The main purpose of their work was to show that the modifications 
of the Mott scattering differential cross section for the scattering of an electron by the Coulomb 
potential of a nucleus in the presence of a strong laser field, can yield interesting physical insights 
concerning the importance and the signatures of the relativistic effects. Their spin dependent 
relativistic description of Mott scattering permits to distinguish between kinematics and spin-orbit 
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coupling effects. Tliey liave compared tlie results of a calculation of the first Born differential cross 
section for the Coulomb scattering of the Dirac-Volkov electrons dressed by a circularly polarized 
laser field to the first Born cross section for the Coulomb scattering of spinless Klein-Gordon particles 
and also to the non relativistic Schrodinger-Volkov treatment. The aim of our work is to provide 
the correct expression for the first-Born differential cross sections corresponding to the Coulomb 
scattering of the Dirac-Volkov electrons. On the one hand, we show that the terms proportional to 
sm(20o) are missing in [T], where 0o is the phase stemming from the expression of the circularly 
polarized electromagnetic field. The claim of [T] that they vanish is not true. These terms do 
not depend on the chosen description of the circular polarization in cartesian components. On the 
other hand, we perform the calculations with some details and throughout this work, we use atomic 
units {h = e = m = 1) where m denotes the electron mass. The abbreviation DCS stands for the 
differential cross section. 

The organization of this paper is as follows : in Section 2, we establish the expression of the 
S-matrix transition amplitude as well as the formal expression of scattering DCS. In Section 3, we 
give a detailed account on the various trace calculations and show that indeed there is a missing 
term proportional to sin(20o) that is not equal to zero. This term as well as a term proportional 
to cos(20o) contribute to (^^^ and multiply the product Js+i{z)Js-i{z), where Js{z) is an ordinary 
Bessel function of argument z and index s. The argument z appearing in the above mentioned 
product will be defined later. Then, we carry out the derivation of the correct expression of the 
scattering DCS associated to the exchange of a given number of laser photons. In section 4, we 
give some estimates of the numerical significance of our corrections. In particular, we compare 
numerically the Dirac-Volkov DCS we have obtained with the corresponding DCS of pp. We end 
by a brief conclusion in Section 5. 

2 The S'-matrix element and the scattering differential cross 
section 

Exact solutions of relativistic wave equations p] are very difficult to obtain. However, in a seminal 
paper, Volkov [3] obtained the formal solution of the Dirac equation for the relativistic electron 
with 4- momentum p'^ inside a classical monochromatic electromagnetic field A^^. These solutions 
are called the relativistic Volkov states. The plane wave electromagnetic field A'^ of 4-momentum 
k^^ {k^k^ = k"^ = 0) depends only on the argument = k.x = k^x^ and therefore A^ is such that 



The 4-vector satisfies the Lorentz gauge condition d^A^^ = or equivalently k^A^ = 0. The 
Dirac-Volkov equation for an electron in an external field is 



where F^^, is the electromagnetic field tensor F^^, = d^A^, — d^A^ and a'^^ = ^[j^, 7*^]. The matrices 
7^^ are the anticommuting Dirac matrices such that 'yt^'j'^ + ^'^^m = 2g^^'^l4^, where g'^'^ is the metric 
tensor g^'^ = diag{l, —1, —1, —1) and I4 is the identity matrix in four dimensions. The solutions of 



A" = A^'ik.x) = A^icP). 



(1) 




(2) 
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Eq.([2]) are the relativistic Dirac-Volkov wave functions 
whith 



and the function S{x) is given by 

S{x) = —p.x 



k.x 



1 



c{k.p) I 2c 



dC (5) 



In Eq.(l3]), u{p, s) represents a Dirac bispinor which satisfies the free Dirac equation and is normahzed 
according to u{p, s)u{p, s) = u*{p, s)'~^^u{p, s) = 2c^. We consider a circularly polarized field 

A = ai cos(0) + 02 sin(0), (6) 

where = k.x. We choose = a"^ = a"^ = A? and 01.02 = 02.01 = 0. The Lorentz condition 
k.A = implies Oi.fc = 02.fc = 0. If one assumes that A^ is quasi-periodic so that its time average 
is zero A^^ = 0, then using the Gordon identity, the averaged 4-current is easily obtained : 

If one sets 
this yields 

Q-Q = Q^Qfj. = mlc^, (9) 



with 



ml = l-^. (10) 



One often calls the averaged 4- momentum a quasi- impulsion. Note that = (Q/c, q). The 
quantity plays the role of an effective mass of the electron inside the electromagnetic field. For 
the study of the process of Mott scattering in presence of a laser field, we use the Dirac-Volkov wave 
functions [3] normalized in the volume V: 

u^) = (11) 



whith 



= 1 + ^^^^ cos(0) + sm(0)) (12) 
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and 

= -^■"-^""^'^) + ^^°^^'^)- 

We turn now to the calculation of the transition amplitude. The interaction of the dressed electrons 
with the central Coulomb field 

^'^ = (-^,0,0,0) (14) 

is considered as a first-order perturbation. This is well justified if Za <^ 1, where Z is the nuclear 
charge of the nucleus considered and a is the fine-structure constant. We evaluate the transition 
matrix element for the transition (i — > /) 



We first consider the quantity 



^0 I 1 .,0 

We have 



i^qf{x)^^^P,^{x) = -^-j==u{p^^ Sf)R{pf)^^R{pMPr. s,)e-^(^(^^'^)-^(^-^«. (16) 



^-i{S{qi,x)-S{q,,x)) ^ exp[i(g/ - qi).x - izsiii{(t) - 0o)], (17) 
where z is such that 



z = ^ai + al (18) 
whereas the quantities cci and 0:2 are given by 

^ {ai-Pi) _ jai-Pf) ^ {a2-Pi) _ (Q2-P/) 
^ c{k.pi) c{k.pf)' ^ c{k.pi) c{k.pf) 

and the phase 0o is such that 0o = arccos(Q;i/2;) — arcsin(a2/;i;) = arctan(Q;2/tti)- It is important 
at this stage to perform intermediate calculations in order to reduce the numbers of 7 matrices that 
will appear when one calculates the scattering DCS. After some algebraic manipulations, one gets 

u{pf,Sf)R{pf)-f^R{pi)u{pi, Si) 
= u{pf,Sf)[Co + Cicos{(j)) + C2sm{(f))]u{pi,Si), (20) 

where the three coefficients Co, Ci and C2 are respectively given by 

Co = -2koa^MPiHPf) 
Ci = c(pi)7°|^^i + c(p/)^i^7° 

C2 = c(p,)7°^^2 + c(p/)^2^7° (21) 
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with c(p) = 2wT^ ko — — cu / c. Therefore, the transition matrix element becomes 

1 



2c{k.p) 



^ ^Z f ,4 _A 

c J ^2Q,V J2QfV 



u{pf, Sf)[Co + Ci cos(0) + C2sm{(f))]u{pi, Si) 



X exp[i{qf — qi).x — izsin{4> — 4>q)]. 
We now invoke the well-known identities involving ordinary Bessel functions Js{z) 



(22) 



1 

cos(0) 
sin(0) _ 



> e 



-izsin(0— </>o) 



-Bis 



> e 



(23) 



with 









hi 







(J,+i(z)e^(*+^)'^° + J,_i(;2)e'(^-^)^°)/2 
(J,+i(^)e^("+^)'^° - J,_i(^)e^(^-^)'^o)/2i ^ 



(24) 



Evaluating the integrals over xq and x yields for 5"/; 



2^ I .1.12 ^^Vi ' 



V2g77J2gyy .ir^o |q/ - - 



(25) 



where the quantity M^f is defined by 



Mjf = u{pf, Sf)[CoBs + CiBu + C2B2s]u{pi, Si). 



(26) 



To evaluate the DCS, we first evaluate the transition probability per particle into final states within 
the range of momentum dq/ 



dWfi 



,ydqf 
"(2^ 



2-72 oo 



E 



T2nS{Qf -Qi- sw) |^(,)|2 V'dq/ 



2QiV.2QfV 
where we have used the rule of replacement 



|q/ - qj - sk|4 



{2ny 



(27) 



[2TrS{Qf - Qi- sw)]'^ 2TrS{0)27rS{Qf - Qi - sw) 

= T27r5{Qf-Qi-sw). 

Next, we have for the transition probability per unit time 

dWfi 



(28) 



dR 



fi 



T 



2v2 oo 



2QiV.2QfV 



E 



27i6{Qf -Qi- sw) ^^{s)^2Vdqf 



|q/ - qi - sk|4 



(27r)3 



(29) 
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Dividing dRji by the flux of incoming particles 



Q^V 



(30) 



then using the relation |q/|(i|q/| 
scattering DCS 

da 



-^QfdQf and integrating over the final energy, we get for the 



dQi 



E 



^ s=oo 

da^'^ 



|q/ - qi - sk| 



Q f=Qi+sw 



E 



dn 



f 



Qf=Qi+sw 



where 



dQi 



Qf—Qi 





q/l 






2 






|q/ - 


- Qi - 


- sk 


4 



(31) 



(32) 



The calculation is now reduced to the computation of traces of 7 matrices. This is routinely done 
using Reduce We consider the unpolarized DCS. Therefore, the various polarization states have 
the same probability and the actually measured DCS is given by summing over the final polarization 
Sf and averaging over the initial polarization Sj. Therefore, the unpolarized DCS is formally given 
by 

da ^ da^' 



dVLi 



E 



d^i 



Qf=Qi+sw 



where 



da^''> 



dn 



f Q f=Qi+sw 

3 Trace calculations. 







1 


C4 


qi 




q/ - qi - sk 
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^eek: 



W|2 



(33) 



(34) 



Qf=Qi 



Since the controversy is very acute and precise about the results of the sum over the polarization 
I J2si J2sf we devote a whole section to the calculations of the various traces that intervene 

in the formal expression of the unpolarized DCS given by Eg. 0341) . We have to calculate 



^EEI<P 



EE ■^/)i'^o5s + CiBu + C2B2s]u{pi, Si) 



with 



^EEl"(p/'^/)^^'^«fe'^* 

Si Sf 



A(^) = [^'^-2koa'MP^)c{Pf)]Bs 



(35) 



(36) 
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Using standard techniques of the 7 matrix algebra, one has 

\ E E \Mf^ P = \Tr{{i>fC + c^)^''\hc + c^)A^^^}, (37) 

Si Sf 



with 



7 iv' "7 
= [7° - 2A:oa2^cfe)c(p;)]5: 

+ [cfe)^2^7° + c(p;)7°^^2]i?2*s- (38) 
There are nine main traces to be calculated. We write them explicitly 

Ml = Tr{{i>fC + c^)C,{i>,c + c^)Co}\Bs\\ 
M2 = Tr{{i>fC + c^)C^{i>,c + ^)Ci}BsBl, 
= Tr{{i)fC + ^)C^{i)iC + ^)C2}B,Bl, 
Ma = Tr{{i>fc + ^)Ci{i)iC + c2)Co}5:Si„ 

M^ = Tr{{i)fC + ^)Ci{i>,c + ^)Ci}\Bu\\ (39) 
^6 = rr{(j^^^c + ^)Ci{i),c + c2)C2}Si.S2*„ 
A<7 = Tr{{i>fc + c2)C2(j^^iC + ^)Co}B2sB:, 
Ms = Tr{{^fc + c')C2{U + c^)Ci}BlB2s, 

M, = Tr{{i)fC + ^)C2{i)iC + c2)C2}|i?2s|'. 

To simplify the notations, we will drop the argument of the various ordinary Bessel functions that 
appear. The diagonal terms give rise to 

A^5 « \Bu? = i(Jf+i + 2Js+iJs-i cos(20o) + Jf-i), (40) 
7U9 oc |S2.|2 = - 2Js+iJs-i cos(20o) + 

So, taking into account the fact that the traces multiplying |-BsP, l-^isP ^^^1 |-B2sP are not zero, 
one expects that terms proportional to J^+i J^-i cos(20o) will be present in the expression of the 
scattering DCS. The first controversy between our work and the result of Szymanowski et al [1] 
concerns the traces Mq andAlg- Since 



M^ (X = |( Jf+i - 2zJ,+i J,_i sin(20o) - 

Ms oc Bl,B2s = f ( + 2«J,+i J,_i sin(20o) - J's-i 



and with little familiarity with the 7 matrix algebra, one can see at once that if the corresponding 
traces are not zero then the net contribution of Mq + Ms will contain a term proportional to 
Js+i Js-i sin(20o)- We shall demonstrate that in what follows. We have 

Me = Tr{{i>fC + ^)Ci{i>,c + ^)C2}BuBl 

= Tr{{^fC + c^) [c(p,)7°^^i + c(p/)^i^7°](Ac + c') 

[c(p.)^2^7° + cipfh'^^2]}B,sBl. (42) 
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From now on, we define a 4-vector 
We can therefore write 



7 



1,0,0,0). 



(43) 
(44) 



Tlien, Eq.(H2!) becomes 

Me 



+ c 



2\}BisBl^. 



(45) 



In [T], tlie autliors claim tliat tfie controversial sin(20o) term disappear because it is proportional to 
terms like Tr{{j[)fC+c^)'~^'^l^(/Li{j[>iC+c^)(^2^1^}- This term as well as Tr{(j(^jc+c^)(^i^7''(jf^jC+c^)7°|if(^2} 
are indeed zero but for Tr{{j[)fC + c^)'^^^(^i{j()iC + c'^)j^lt^2} and Tr{(/^/c + c^)^ilt'y^{^iC + c^)(/i2^'^^} 
this is no longer true. These terms are not zero and we give explicitly their values 

Tr{{i>fC + c')^'Ui{i>iC + c')^'m = Tr{{i>jc + c')i^^^\i>,c + c^)<li2H] 

= Sw'^{{ai.pf){a2.pi) + {ai.pi){a2.pf)}. (46) 

In most cases, the various traces are zero except when the cyclic process of taking scalar products 
of pairs comes to products such that 



{k.r]){k.r]){ai.pi){a2.pf), 
{k.7]){k.r]){ai.pf){a2.Pi), 



(47) 



in which case, one has contributions proportional to w'^{ai.pi){a2-Pf) and w'^{ai.pf){a2-Pi) respec- 
tively. Explicitly, we give the result for A4e and Ais- 
One has 



Mf 



w 



{2 sin(20o 



ai.pi) {a2.pf) ^ {a2-Pi) (ai-Pf) 



(k.pi) (k.pf) (k.pi) {k.pf 
+i[-{{ai.pi){a2.pf) + {ai.pf){a2.pi)}J^_^ 
+{{ai.pi){a2.pf) + {ai.pf){a2.pi)}Js+i]}, 



Js+lJs-l 



(48) 



while A^8 is given by 
-Ms 



w 

— {2sin(20o 



(ai.pi) {a2.pf) _^ {a2-Pi) (oi-P/) 



(k.pi) (k.pf) {k.pi) {k.pf 
-i[-{{ai.pi){a2.pf) + {ai.pf){a2.Pi)}J^_i 
+{{ai.pi){a2.pf) + {ai.pf){a2.pi)}J]+i]}. 



Js+lJs-l 



(49) 



The fact that complex numbers appear in the expressions of and Aig is not surprising since 
the former is the complex conjugate of the latter and their real sum is such that 



4w^ 



sin(20o) 



[ai-Pi) {O'2-Pf) _^ {a2-Pi) (ai-Pf) 



_ (k.pi) (k.pf) (k.pi) (k.pf) _ 



Js+lJi 



(50) 
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So, the first controversy is settled and there is indeed a term containing sin(20o) in the expression 
of the scattering cross section. We have written a Reduce program that calculates analytically the 
traces in Eq. fl37p . Before writing our Reduce program, we have extensively studied the textbook by 
A. G. Grozin [5j which is full of worked examples in various fields of physics particularly in QED. 
We give the final result for the unpolarized DCS for the Mott scattering of a Dirac-Volkov electron 



2 







1 1 


C2 






|q/ - qi - sk| 


4 



{J^ji + U'Li + Jli)B + J,_,)C + J,( + J.+i)£>} , (51) 



where for notational simplicity we have dropped the argument z in the various ordinary Bessel 
functions. The coefficients A, B, C and D are respectively given by 



c^{k.qf){k.qi) \{k.qi) {k.qj) ^ 

^ {a^fu'^ / {ai.qf) {ai.qi) {a2-qf) {a2-qi) \ 

2c^{k.qf){k.q,) 2c^\{k.qf) {k.q,) {k.qj) (k.q,) J 2 



a' Ak.qf) jk.q,) aW « ^ _ , 

A\k.q,) + {k.qjY 2c2(M/)(A:.gO^^ ^ 

''"(^/-Q^)f77^-77^V (53) 



2c^ \{k.qf) {k.q. 



C = -^^^j—^^j—^{cos(2<i)Q){{ai.qf){ai.qi) - (aa.g/) (aa.gi)} + 

sin(20o){(ai.g/)(a2-gi) + (ai.gi)(a2.g/)}), (54) 

u ( Q^{A.qf) , Qf{A.q.i) \ 

c {-jk:^ + iMir j ' ^ ^ 

where A = ai cos(0o) + a2 sin(0o). 

3.1 Comparison of the coefficients. 

The argument about the missing term proportional to sin(20o) having been given a convincing 
explanation, we now turn to other remarks along the same lines since there are indeed other dif- 
ferences between our result and the result of [I]. We discuss now the difference occurring in our 
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expression of the coefficient A and the corresponding one of p^. In their expression multiplying 
the product 2J^(^), the single term should come with a coefficient |. We have written a 

second Reduce program that allows the comparison between the coefficient A of [Ij and the coeffi- 
cient A of this work. There are so many differences between our result and the result they found 
for the coefficient B that we refer the reader to our main Reduce program [8]. The coefficient C 
has already been discussed. As for the coefficient D, we have found an expression that is linear in 
the electromagnetic potential. In a third Reduce program, it is shown explicitly that if we ignore 
the first term in the coefficient multiplying Js{Js-i + Js+i) given in [1], one easily gets the result we 
have obtained. This term does not come from the passage from the variables {p, p) to the variable 
{q,q)- The introduction of such 4- vector q is not useful, makes the calculations rather lengthy and 
gives rise to complicated expressions. As a supplementary consistency check of our procedure used 
in writing the main Reduce program, we have reproduced the result of the DCS corresponding 
to the Compton scattering in an intense electromagnetic field given by Berestetzkii, Lifshitz and 
Pitaevskii [6]. 



4 Results and discussion 
4.1 Kinematics of the collision 

For the description of the scattering geometry, we work in a coordinate system in which k| je^. This 
means that the direction of the laser propagation is along the Oz axis. To avoid any confusion, we 
will compare the Dirac-Volkov DCS (26) of pLj with the corresponding DCS (IMl) we have obtained in 
the same coordinate system. The spinless DCS will also be discussed as well as the non relativistic 
one. We begin by defining our scattering geometry. In our system, the vector Pj is such that [e^;, 
meaning that the undressed angular coordinate of the incoming electron are 6i = 90°, 0j = 0°. For 
the scattered electron, the vector pj is such p^ G (ey,e^), meaning that —180° < 9f < 180° and 
(pf = 90°. With this choice, we have been able to reproduce qualitatively the results and all the 
figures of [1]. The reason underlying this choice of the coordinates is the following. The angles 
{6i, (pi) of Pj (the same holds for p^) are the intrinsic angular coordinates of the incoming electron. 
As the vector is defined through p^ via the relation 

= - w^^^ 

we cannot define intrinsic angular coordinates using qj. When the electron is subjected to the radia- 
tion field, it acquires new angular coordinates that can easily be determined. The key quantity that 
gives an idea of the dependence of qj(and qj) on the spatial orientation of the electron momentum 
due to (fc.pj)(and (k.pf)) is the cosine of the angle between q^ and qj. While 

Pi-P/ = |Pi||P/l cos(pI7p/), (57) 

with 

cos(pJ7p/) = sin{9i)sin{6 f)cos{(j)i — (j)f) + cos{9i)cos{9f), (58) 

we have 

Qi-q/ = |qi||q/l cos(q^/), (59) 
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with 



I II I / — A I II I / — N a'^w f\pi\cos{ei) |p/| cos(6'/)\ 
|q,||q,|cos(q.,q,) = |p.| |p,| COs(p., p,) - ^ + ^^^^^^ j 

[a_j_w 

A{k.pf){k.p,)c^ ^ ^ 

From these relations, one deduce that in the hmit of low incoming electron energies and moderate 
field strength, cos(q7rqf) and cos(p7rPf) are very close therfore the second and third term of the 
RHS of Eg. (1601) can be safely neglected. For high incoming electron energies and intense field 
strength, the difference between the two cosines increases and these terms cannot be neglected. We 
shall give for the sake of illustration, tables that compare these two cosines for the three regimes we 
shall investigate, namely the non relativistic-moderate field strength regime, the relativistic-strong 
field strength regime and finally the relativistic-intense field strength regime. We choose the same 
value as pQ for the laser angular frequency w = 0.0430 a. m for all the numerical calculations. This 
typical near infra-red angular frequency is that of a neodymium laser. The other parameters are 
the electric field strength e and the relativistic parameter 7 = — = 1/ \Jl — jc?. This 

parameter fixes the incoming electron total energy Ei via the relation Ei = 0^7 from which one 
deduces the corresponding incoming electron kinetic Tj by subtracting the rest energy (c^ in a.u) 
Ti = Ei — = c^Ij — 1). Before beginning our discussion, we would like to make general comments 
on the figures obtained in [Ij starting with Figure [3l This figure does not represent the envelope of 
the controversial generalized equation (26) of that work. Indeed, we shall see that it represents the 
envelope of the non relativistic DCS given by Eq.(34) in [T]. We give the correct envelope for the 
relativistic calculations obtained by using either Eq. (!M|) of our work or Eq.(26) of [Ij. In Figure ([6]) 
of [1], there is a difference between the Dirac-Volkov DCS (26) and the spinless particle DCS (30) 
though the overall behaviour is smoothly oscillatory. The results we have obtained show the same 
oscillatory behaviour. The curves for the Dirac-Volkov DCS (26) of [1] and the Dirac-Volkov DCS 
( IMl) of our work are almost identical while the difference between two relativistic DCSs and the 
spinless particle DCS given by Eq.(30) of [Ij is less important than in Figure 6 of [Ij. Figure 7 of is 
the only figure we agree with. In Figure 8 of we disagree with the behaviour of the Dirac-Volkov 
DCS (26) of [Ij particularly for small angles around 6f = 0°. When programing Eq.(26) of [1], we 
obtained a value for the Dirac-Volkov DCS at ^/ = 0° of nearly 3.2 10~^^ a.u instead of the 2.2 10"^'* 
a.u indicated in Figure 8 of [T]. Moreover, the electric field strength e being a key parameter (as 
well as the incoming electron total energy), we have compared our Dirac-Volkov DCS and the Dirac- 
Volkov DCS (26) of [1] and we have come to the following important conclusions. First, for the 
non relativistic and low and low-intensity field strength regime (7 = 1.0053 a. m, e = 0.05 a. m) and 
for the relativistic regime and increasing field strength (7 = 2.00 a. m, e = 1.00 a. m) the differences 
between our results and the results found in [T] are small but approach one percent. Second, we 
have a different picture for the relativistic-high intensity regime (7 = 2.00 a. m, e = 5.89 a.u) where 
the missing terms in lead to values of the Dirac-Volkov DCS (26) of [1] that over-estimate the 
corresponding DCS flM|) of our work. Even in the non relativistic regime (7 = 1.0053 a.u) but for 
increasing field strength, the difference between our results and the results of [1] begins to appear 
clearly. 

We turn now to a qualitative and quantitative discussion of the physical process. We shall comment 
and analyze the results obtained in [Ij in the light of those we have obtained bearing in mind that 
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we can hardly escape rephrasing the physical insights and explanations contained in [T]. Our 
disagreement is quantitative since we have shown in the first part of this work that the expression 
(26) of [T] contains errors and a missing term proportional to sin(20o)- So, our primary task is 
to assess the importance of this errors and missing term and to what extent they modify the 
quantitative and qualitative contents of 

4.2 The non relativistic-low electric field strength regime 

In this regime, we choose as in pj 7 = 1.0053 a.u for the relativistic parameter and e = 0.05 a.u for 
the electric field strength. This relativistic parameter corresponds to an incoming electron kinetic 
energy Tj = 100 a.u = 2.700 kev. With our choice of the angular parameters, we compare in Table. 
[Usome values of cos(pi7p/) and cos(q^7q/)- 



cos(prP/) 


cos(q;7q/) 


0.853553 


0.853631 


0.850868 


0.850942 


0.847923 


0.847993 


0.844719 


0.844786 


0.841259 


0.841322 



Table 1: 

The difference is small so in this regime and the coordinates of the undressed electron are nearly the 
same as that of the dressed electron. We plot in the upper part (a) of Figure [1] the non relativistic 
DCS given by Eq.(34) of [T] and in the lower part (b) of the same figure, the generalized Dirac- 
Volkov DCS given either by Eq.(26) of [Ij or Eq.( |Ml) of our work as a function of the final electron 
energy scaled to the photon energy. The scattering angle is large enough so that an important 
number of photons can be exchanged in the course of the collision. In this low-intensity regime, the 
envelope of the non relativistic DCS is qualitatively different from the envelope for the Dirac-Volkov 
and Klein-Gordon DCSs. 

The observed cutoffs occur at Smin = —522 and Smax = 582 for the non relativistic DCS and 
Smin = —474 and Smax = 474 both for the Dirac-Volkov and Klein-Gordon DCSs since the argument 
that appears in the ordinary Bessel functions is the same for both DCSs. So the comments made 
in [1] concerning the interpretation of the envelope obtained do not apply for the Dirac-Volkov and 
Klein-Gordon cases. While the spectrum of Figure (l.a) of our work (which is identical to that of 
Figure (l.a) of |1J) exhibits an overall asymmetric envelope with peaks of negative energy transfer 
higher than peaks of positive energy transfer, this asymmetry is less pronounced in the case of the 
Dirac-Volkov and spinless particle DCSs. This emphasized asymmetry in the non relativistic case 
can easily be traced back by a close look at Eq.(34) of [Ij. Indeed, the non relativistic DCS depends 
on Jg{z) {z depends only weakly on s) so the asymmetry can only come from the dependence of 
the modulus of the final momentum qj on the number of the transferred photons s according to 
Eq.(25) of yy. We explicitly write this equation (with our notation) 

II A |2 , o ^'^i , 2^^V^^ ff,-i\ 

|q/| = Mqil +2s— - + s — . (61) 
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Figure 1: (a): Envelope of the non relativistic differential cross-section da/dQ scaled in unit of 10~^ 
a.u as a function of energy transfer Qf — Qi scaled in units of the laser photon energy w for an 
electric field strength of 0.05 a. m. and a relativistic parameter 7 = 1.0053 a. m, (b). Envelope of the 
relativistic differential cross-section {da/dVL)i:,v scaled in unit of 10"'' a.u as a function of the laser 

photon energy w for the same parameters. The envelope for ((icr/(if2)§|/ and {da / dVL)KG are almost 
identical 

Also, the denominator Iq,; — qj + sk|^ depends linearly on the number s of transferred photons 
and this gives rise to an asymmetric envelope for the non relativistic DCS. In Figure 2, we plot 
the behaviour of |qi — q/ + sk|~^ as a function of the number s of photons transferred. As 




-100 -50 50 100 



Pfioton energy transfer 

Figure 2: Behaviour of |qj — q/ + sk|~'^ as a function of the number s of the photons transferred 
with 7 = 1.0053 and e = 0.05 a.u 

mentioned in [1], we have an enhancement of negative over positive-energy transfer cross-section. 
The DCSs fall of abruptly beyond the points where the argument of the Bessel functions equal to the 
order. For the Dirac-Volkov and Klein-Gordon DCSs, this cutoff occurs ( up to machine precision) 
numerically for s = ib474 and an argument z of the ordinary Bessel functions almost constant and 
equal to 380.016. However Figure (l.b) shows a visual cutoff for s = ±392 since the infinitesimal 
contributions to the DCSs cannot be plotted. For the non relativistic DCS, the numerical cutoff 
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occurs (again up to machine precision) for s = —606 and s = 685. The visual cutoff occurs for 
s = —500 and s = 590. The difference As = 20 between our resuhs and that of [I] is just a matter 
of convention. We now analyze the angular distributions. We have summed as in ±100 peaks 
around the elastic one in order to draw the angular dependence of the DCS. In Figure 6. of [1], 
the accumulated DCS is shown for an electric field strength e = 0.05 a. u. The computer code we 
have written calculates the Dirac-Volkov DCS flMl) of our work, the Dirac-Volkov DCS (26) of [1], 
the spinless particle DCS and the non relativistic DCS. At least, in the non relativistic regime, our 
results and that of [1] agree very well and are both close to the results for a spinless particle. We 
give in Figure [3] the angular distribution of the various DCSs. 




angle 9^ (degree) 



Figure 3: Summed differential cross sections da/dQ of ±100 peaks around the elastic one as a 
function of the angle Of for a relativistic parameter 7 = 1.0053 a. u and an electric field strength 
e = 0.05 a.u. The solid line denotes the result for Dirac-Volkov electrons, the long dashed one 
sketches the values for {da/dQ)^^Y and the short dashed is the result for spinless particles 

Apart from minor differences, all three calculations exhibit maxima for 6f = 0° and ±180°, a 
giggling oscillatory behaviour (as in [1]) and minima slightly shifted from ±90° ( at ±84°). Let aside 
the order of magnitude, we have in our case, three DCSs that are close to each other and not as 
differentiated as shown in Figure 6. of [1]. We do not agree at all with the results shown in Figure 
6. In particular, our extrema for the various DCSs are shown in table ([2]) (scaled in 10~^ a.u). 





{da/dil)Dv 


(da/dny^^y 


{da/dn)KG 


-180° 


4.32027 


4.32027 


4.34311 


-84° 


3.01486 


3.01486 


3.03079 


0° 


4.326 


4.326 


4.34886 


84° 


3.01486 


3.01486 


3.03079 


180° 


4.32027 


4.32027 


4.34311 



Table 2: 
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So, this adds to the controversy. Even if we use the expression for the Dirac-Volkov DCS given 
by Eq.(26) of p], we have a different figure for the non relativistic regime. If we now increase the 
electric field strength from e = 0.05 a. m to e = 1.00 a. m, the agreement remains good between the 
three relativistic calculations. There is still a maximum at 6'/ = 0° while the minima are shifted 
towards ±117°. To give an idea the small differences between our result and the result of [1] for the 
Dirac-Volkov DCS, we have plotted in the upper part (a) of Figure 4, the ratio of the DCS given by 
Eq.(26) of ppj to the DCS given by Eq. flMl) of our work as a function of the angle Of for e = 1.00 u.a. 
The ratio R is defined by 

/ da 



R 



Idv 



' _da_ 



(62) 



}DV 



The deviations from the expected value 1 are shown and have the same shape as the corresponding 




angle 6^ (degree) 



(b) 



90 135 
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Figure 4: (a) Ratio R of the two Dirac-Volkov DCS for 7 = 1.0053, e = 1.00 a. u and s 
(b): Ratio R of the two Dirac-Volkov DCS for 7 = 1.0053 a. m, e = 5.00 a. u and s = ±100 



±100. 



DCS. However, for increasing electric field strength, the values for this ratio are not close to 1. For 
a relativistic parameter 7 = 1.0053 a. u and for an electric field strength e = 5.00 m. a and s = ±100, 
our results for the Dirac-Volkov DCS and the corresponding results of [1] do not agree at all. In 
the lower part (b) of Figure (4), there is an over estimation varying from 22.5% to 30% with some 
peaks giving an over estimation of up to 45% for the DCS (26) of pLj compared to the corresponding 
DCS ( 151|) of this work. All these peaks are nearly multiples or submultiples of an angle close to 
7r/4. 



4.3 Relativistic-strong electric field strength regime 

For the relativistic regime, we have chosen the parameters of [1] 7 = 2 which corresponds to an 
incoming electron total energy Ei = 2c^ or a Tj = 0.5116 Me\^. The electric field strength is now 
e = 1.00 a. M. Some cosines of the angles cos(qi7q/) and cos(p^7p/) are shown in Table ([3]) 
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cos(p^/) 


cos(q^/) 


0.853553 


0.855757 


0.850868 


0.852775 


0.847923 


0.849543 


0.844719 


0.846062 


0.841259 


0.842332 



Table 3: 

In this regime, dressing effects are important. The envelope of the energy distribution of the 
scattered electrons is similar to the one displayed in the lower part (b) of Figure [H However, there 
is a more important asymmetry than in the non relativistic regime with was to be expected. The 
corresponding cutoffs are (-170000, 194000) for ( Jf^) and (-66000, 64000) for ( J^) ,( J^)™ 
and ( J^)^^- The three relativistic calculations lead to angular distributions peaked in the direction 
of the laser propagation Of = 0°. The two Dirac-Volkov DCSs (solid line and long dashed line) are 
slightly different only in the vicinity of the two minima located at Of ~ ±33°. In this regime, 
the shape of the ratio R is similar to that of the corresponding DCSs. This ratio is equal to 1 
for 6f = 180° but there is now an overall amplitude of 8.10"^ around the expected value 1. If we 
increase the electric field strength from e = 1.00 a.u to e = 5.89 a.u and keep the same value of the 
relativistic parameter 7 = 2, the difference between our Dirac-Volkov results and the corresponding 
results of [1] becomes important. 

In the upper part (a) of Figure 5, we give the two Dirac-Volkov DCSs and in the lower part (b) of 
the same figure we give the ratio R of the two DCSs for the relativistic parameter 7 = 2.00 and an 
electric field strength e = 5.89 a.u. For the angles 6f = ±180°, the ratio is i? ^ 1.01 while for the 
peak in the direction of the laser propagation, 6f = 0°, the ratio is i? ~ 2.34. 

4.4 Relativistic and high electric field strength regime 

To study this regime, we use the same parameters as in [Ij, that is a relativistic parameter 7 = 5 
or an incoming electron kinetic energy Tj = Ac^ a.u = 2.045 Mev. The electric field strength is 
e = 5.89 a.u. In table(jl]), some cosines of the angles cos(qi7qf) and cos(pi7Pf) are given. In this 
regime, the dressing of the electron angular coordinates is very important. 



cos(p^/) 


cos(qrq/) 


0.853553 


0.866661 


0.850868 


0.861677 


0.847923 


0.856549 


0.844719 


0.851270 


0.841259 


0.845833 



Table 4: 
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Figure 5: (a) Summed differential cross sections da/dQ of ±100 peaks around the elastic one 
as a function of the angle 6f for a relativistic parameter 7 = 2.00 and an electric field strength 
e = 1.00 a. M. The solid line denotes the result for the Dirac-Volkov electrons, the long dashed one 
sketches the values for {da/dQ)^^y and the short dashed is the result for spinless particles, (b): 
Ratio R of the two Dirac-Volkov DCSs for a relativistic parameter 7 = 2.00 a. m and an electric field 
strength e = 5.89 a. u 

In the upper part (a) of Figure 6, we show the various DCSs. For angles 9f = ±180°, the agreement 
between our results and the results of [T] is good but deteriorates for small values of For 9f = 0°, 
the result of our work gives a value ( scaled in 10"^"^ a. m ) {■^)dv = 1-15 while de corresponding 

result found using Eq.(26) of [Ij is ( — 3-204 . Our results (solid line) are always smaller than 
the results for spinless particles while those obtained using Eq.(26) of P are greater than {-^)kg 
for small angles around the direction of the laser propagation. In the lower part (b) of Figure 6, we 
show the ratio R defined by Eq. (!62|) . For 9f = 0°, this ratio is R = 2.787. 

5 Conclusion 

In this work, we derived the correct expression of the first Born differential cross section for the 
scattering of the Dirac-Volkov electron by a Coulomb potential of a nucleus in the presence of a 
strong laser field. We have given the correct relativistic generalization of the Bunkin and Fedorov 
treatment ^ that is valid for an arbitrary geometry. We are adamant that the core of the whole 
controversy stems from the fact that in pQ, the vector r]'^ introduced in Eq. fl43l) of our work has not 
been properly dealt with while it is the common method to use when a trace contains a 7° matrix. 
Any standard QED textbook introduces this very elementary method. Comparison of our numerical 
calculations [9J with those of Szymanowski et al. [1] shows qualitative and quantitative differences 
when the incoming total electron energy and the electric field strength are increased particularly 
in the direction of the laser propagation. The difference between our results and those of jl] can 
only be traced back to the mistakes and the omitted term in Eq.(26) of fTj. The corrections that 
we made allowed us to study other processes that were published in Physical Review A., namely an 
first article concerning the relativistic electronic dressing in laser assisted electron-hydrogen elastic 
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Figure 6: (a): Summed differential cross section da/dQ of ±100 peaks around the elastic one as 
a function of the angle 6f for a relativistic parameter 7 = 5.00 a. m and an elastic field strength 
e = 5.89 a.u. Th solid line denotes the result for Dirac-Volkov electrons, the long dashed one 

sketches the values for {da/dQ)'i)y and the short dashed is the result for spinless particles, (b): 
Ratio R of the two Dirac-Volkov DCSs for the same values of the relevant parameters, 7 = 5.00 a.u 
and e = 5.89 a.u 



collisions [TU], another concerning the process of Mott scattering in an elliptically polarized laser 
field [TT] as well as a third work dealing with the process of Mott scattering of polarized electrons 
in a strong laser field [12]. For the difficult process of ionization of atomic hydrogen by electron 
impact, we published an article concerning the importance of the relativistic electronic dressing in 
laser-assisted ionization of atomic hydrogen by electron impact [T3]. All these works relied heavily 
on the corrections that we made in this work. 
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